function Idlist_compose =  mySSA_chooseRCs(Lambda,RCs,Fs,threshold)
    % Idlist_compose =  mySSA_chooseRCs(Lambda,RCs,Fs,threshold)
    %
    % 说明：判断mySSA_decompose分解的重建因子分别的主要频率和解释率，返回重建
    %       因子的索引矩阵（逻辑矩阵），进行重建。
    % 输入：
    %       Lambda： double列向量，特征值
    %       RCs：    double列向量构成的矩阵，每个列向量是一个重建因子
    %       Fs：     采样频率
    %       threshold： 选用的重建因子要在什么样的频率范围内
    % 输出：
    %       Idlist_compose：向量，索引矩阵，选定的用于重建的因子。
    
    [Idlist_95,Idlist_99] = sub_plotContribRate(Lambda);

    [RCs_Fs,RCs_P] = sub_getDominantFrequency(RCs,Fs,Idlist_99);
    sub_plotDominantFrequency(RCs_Fs,RCs_P,Lambda,Idlist_95,Idlist_99);
    Idlist_Fs = sub_chooseDominantFrequency(RCs_Fs,threshold);
    Idlist_compose = Idlist_95 & Idlist_Fs;
end